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Abstract. A particle-number reprojection method is applied in the framework of the 
shell model Monte Carlo approach to calculate level densities for a family of nuclei 
using Monte Carlo sampling for a single nucleus. In particular we can also calculate 
J> \ level densities of odd-even and odd-odd nuclei despite a new sign problem introduced 

by the projection on an odd number of particles. The method is applied to level 
densities in the iron region using the complete pf + f/g/ 2 -shell. The single-particle 
level density parameter a and the backshift parameter A are extracted by fitting the 
microscopically calculated level densities to the backshifted Bethe formula. We find 
good agreement with experimental level densities with no adjustable parameters in the 
microscopic calculations. The parameter a is found to vary smoothly with mass and 
does not show odd-even effects. The calculated backshift parameter A displays an odd- 
even staggering effect versus mass and is in better agreement with the experimental 
data than arc the empirical values. 
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INTRODUCTION 

The nucleosynthesis of heavy elements takes place by radiative capture of neu- 
trons (s and r process) and protons (rp process) in competition with beta decay. 
In the statistical regime, neutron and proton capture rates are proportional to the 
level density of the corresponding compound nucleus [1] . Most theoretical models of 
level densities are based on the Fermi gas model, e.g., the Bethe formula [2], which 
describes the exponential increase of the many-particle level density with both ex- 
citation energy and mass number. In the backshifted Bethe Formula (BBF) (see, 
e.g. Ref. [3]), shell corrections and two-body correlations are taken into account 
empirically by introducing a backshift A of the ground state energy. The BBF of- 
fers a good description of the experimentally determined level densities when both 
a and A are fitted for each nucleus [4,5]. The overall systematics of a and A were 
studied empirically but it is difficult to accurately predict these parameters for a 
particular nucleus. 



The interacting shell model takes into account both shell effects and residual 
interactions and constitutes an attractive framework for calculating accurate level 
densities. However, conventional diagonalization methods are limited by the size 
of the model space. Full major shell calculations are presently restricted to nuclei 
with A < 50 (in the p/-shell) [6] . The development of quantum shell model Monte 
Carlo (SMMC) methods [7,8] allows the calculation of finite and zero-temperature 
observables in model spaces orders of magnitude larger than those that can be 
treated by conventional diagonalization techniques. Recently the SMMC method 
was successfully adapted to the microscopic calculations of nuclear level densities 



The applications of fermionic Monte Carlo methods are often limited by the so- 
called sign problem, which causes a breakdown of the method at low temperatures. 
A practical solution was developed in the nuclear case [8] , but the resulting extrap- 
olation errors were found to be too large for accurate calculations of level densities. 
Instead we have constructed good-sign interactions that include the dominating 
collective components of effective nuclear interactions [10] and were proven to be 
realistic for the calculation of level densities. 

The SMMC method is based on a representation of the many-body imaginary- 
time propagator as a functional integral over one-body propagators in fluctuat- 
ing auxiliary fields, known as the Hubbard-Stratonovich transformation [11]. The 
many-dimensional integration is then performed by Monte Carlo. The SMMC 
method is computationally intensive. In particular, level density calculations re- 
quire computation of the thermal energy at all temperatures. If this procedure is 
to be repeated for a series of nuclei, the calculations quickly become very time- 
consuming. Recently we introduced a novel particle-number reprojection method 
[12] with which we can calculate nuclear observables for a series of nuclei using 
the Monte Carlo sampling for a single nucleus. The weight function used in the 
sampling is proportional to the partition function of a fixed even-even or N = Z 
nucleus. Thermal observables for neighboring nuclei are then calculated by repro- 
jection on different particle numbers (both even and odd). This technique offers 
an economical way of calculating level densities for a large number of nuclei. 



THE SHELL MODEL MONTE CARLO METHODS 



A general many-body Hamiltonian containing up to two-body interactions can 
be written in the following quadratic form: 



where p a are one-body densities. Using the Hubbard-Stratonovich transformation, 
the imaginary-time many-body propagator e~^ H can be represented as 



[9] 
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where G(a) is a Gaussian weight and U a is a one-body propagator of non-interacting 
nucleons moving in fluctuating time-dependent auxiliary fields er(r). The canonical 
expectation value of an observable O at inverse temperature (3 is calculated from 

~ J D[v]G(v)Tr A (dU a ) 
{U)A - J D[a]G(a)Tr A U a ' {6) 

where Tr_4 is a canonical trace in the subspace of fixed particle number A. 
In practice we project on both neutron and proton number, N and Z, re- 
spectively, so A denotes the pair (N,Z). Introducing the notation (X a ) w = 
J D[a]W(a)X a / J D[a]W(a), where W(a) = G(a)Tr A U a , Eq. (3) can be written 

as 

(O)a = (Tr A (OU a )/Tr A U a ) w . (4) 

For a good-sign interaction and for even-even or N = Z nuclei, the weight function 
W(a) is positive-definite. The a-fields are sampled according to W(a) and thermal 
observables are calculated from (4). 



THE PARTICLE-NUMBER REPROJECTION METHOD 

We assume that the Monte Carlo sampling is done for a nucleus with particle 
number A. The ratio Z A i/Z A between the partition function of another nucleus 
with particle number A' and that of the original nucleus A is written as 

Z A ,(P) _ r Fr A ,e-P H /Tr A ,U ax 

(5) 



Z A (/3) Tr A e~P H \Tr A U a / w 

The expectation value of an observable O for nucleus with A' particles is calculated 
from 

/ / Tr^duA / Tr A rUcr \ 

\ Q I* ~ / Tr 4 ,l/ CT \ • ( 6 ) 

\ TtaU* /w 

The Monte Carlo sampling is carried out using the weight function W{a) which is 
proportional to the partition function of nucleus A, and Eq. (6) is used to calculate 
the same observable for nuclei with A 1 ^ A. 

In the calculations of level densities we used the Hamiltonian [9] 

H = E e a n a + goPlW ■ - x £ k x O^ ■ O™, (7) 

a A=2 



where 



^ (A ' T)t = ^nja\\Y,\\j b )[alxa] b }^\ 

ab 
ab 

The modified annihilation operator is defined by dj :rn:fnt = (— l) J_m+ 5~ mt a 7) _ mi _ mt , 
and a similar definition is used for p( A ' T ). The single-particle energies e a are cal- 
culated in a central Woods-Saxon potential V(r) plus spin-orbit interaction. g 
is a monopole pairing interaction strength determined from experimental odd- 
even mass differences. The quadrupole, octupole and hexadecupole interaction 
terms in (7) are obtained by expanding a separable surface-peaked interaction 
v(r,r') = —x(dV/dr)(dV/dr')5(r — r') [13] whose strength x is determined self- 
consistently. The parameters k 2 = 2, fc 3 = 1.5 and k 4 = 1 are renormalization 
constants that take into account core polarization effects. Both the pairing and the 
surface-peaked interactions are attractive and therefore have a good sign [8] . 

In the particle-number reprojection method described above we have assumed 
that the Hamiltonian H is independent of A. Suitable corrections should be made 
if some of the Hamiltonian parameters vary with A. In the iron region we find 
that x depends only weakly on the mass number A, and the pairing strength go 
is constant through the shell. The largest variation in this mass region is that of 
the single-particle energies. The thermal energy of a nucleus with A' = (N 1 , Z') 
can then be estimated from E A '(fi) ps J2 a [ e a(A') — € a (A)]{n a )^ + (H) A r, where H 
is the Hamiltonian for a nucleus with A = (N, Z) . This estimate for E^i is an 
approximation since we are still using the propagator e~@ H with the Hamiltonian 
H for nucleus A (instead of A'). This is a good approximation if we reproject on 
nuclei with N' — Z' values close to N — Z (the Woods-Saxon potential depends 
on iV — Z). In the applications below this leads to negligible errors in the level 
densities. 



(8) 



APPLICATIONS AND RESULTS 

In this section we present applications of the particle-number reprojection 
method to nuclei in the iron region. Since we are interested in level densities around 
the neutron and proton resonance energies we use the complete (pf + pg^-shell. 
This model space contains both positive and negative parity states. 

We perform the direct Monte Carlo sampling for the even-even nucleus 56 Fe and 
the N = Z nucleus 54 Co (both have a good sign for the interaction (7)). The 
thermal energies of 53_56 Mn, 54_58 Fe and 54_60 Co were reprojected from 56 Fe, while 
those of 50 ~ 52 Mn and 52 > 53 Fe were reprojected from 54 Co. The calculations were 
done for (3 values up to 2.5 MeV^ 1 . At small (5 (< 1) the calculations were done 
in a smaller step of Af3 = 1/16. At larger (3 points, the Monte Carlo calculations 
become more time-consuming, and we doubled our step size to A/3 = 1/8. For 
each (3 point we took about 4000 independent samples. The reprojected energies 
usually have larger statistical error at large values of (3. To calculate reliable ground 
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FIGURE 1. The SMMC thermal energies vs. (3 for 50 ~ 56 Mn. Shown on the far right are the 
extrapolated ground-state energies. The inset shows the SMMC thermal energies (diamonds) for 
54 Mn at large [3 values. The triangles are the energies obtained by averaging the large-/3 results 
above the corresponding [3. The extrapolated ground state energy is shown by a horizontal line. 



state energies we performed direct Monte Carlo runs for some of the reprojected 
nuclei at several values of (3 between 1.75 and 2.5. For (3 > 2.5 the statistical 
error for the thermal energy of an odd-even nucleus becomes too large to be useful. 
Since the thermal energy of an odd-even nucleus is already close to its asymptotic 
value at these large f3 values, we could extract the ground state energy to within 
an accuracy of ~ 0.3 MeV. 

Fig. 1 shows the SMMC thermal energies versus (3 for manganese isotopes. The 
staggering observed in the spacings of the thermal energies at large (3 is a pairing 
effect. The inset of Fig. 1 shows the SMMC thermal energies of 54 Mn (diamonds 
with error bars) at large values of f3. It demonstrates the procedure we used to 
extract the ground state energy. 

The level density is related to the partition function by an inverse Laplace trans- 
form 

PAE) = r yLeP E *Z A .{P). (9) 
J-ioo 2m 

The partition function Zj^{(3) is computed from the SMMC thermal energies by 
integrating the thermodynamic relation — d In Z^/d/3 = E^/(j3). The average level 
density is then calculated by evaluating (9) in the saddle point method 

p* = {2it(3- 2 C A ,y l/2 e s * , (10) 

where Sa' = (3Ea> + lnZ^/ and Cai = —f3 2 dEA'/df3 are the canonical entropy and 
heat capacity, respectively 
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FIGURE 2. The SMMC level densities of the 50 " 56 Mn isotopes. The solid lines are the fit of 
the calculated level densities to the BBF (11). 

Fig. 2 shows the level densities for the manganese isotopes of Fig. 1 as a function 
of excitation energy. These densities are fitted to a modified version of the BBF [3] 



Here t is a thermodynamic temperature determined by the relation E x — A = at 2 —t. 
Eq. (11) differs from the usual BBF in the term t which appears in the pre- 
exponential factor, and gives a better fit to the SMMC level densities at lower 
excitation energies. The solid lines in Fig. 2 are the fitted BBF level densities of 
Eq. (11). The fitting is done in the energy range E x < 20 MeV and is usually good 
down to ~ 1 MeV for even-even nuclei (for which A is positive), or even below 
1 MeV for odd- A nuclei. The reduced pairing correlations in odd-odd nuclei are 
clearly observed in the level densities of Fig. 2. The backshift parameter A for the 
odd-odd nucleus 54 Co is lower than A for the odd-even nucleus 55 Co, leading to a 
higher level density for 54 Co despite its smaller mass. 

The level density parameters a and A were extracted by fitting Eq. (11) to the 
microscopic SMMC level densities, and are shown in Fig. 3 versus mass number 
A. The SMMC results (solid squares) are compared with the experimental data 
(x's) quoted in Refs. [4] and [5]. The solid lines describe the empirical formulae 
of Refs. [14]. The SMMC values of a depend smoothly on the mass A, unlike 
the values predicted by the empirical formulae. The pairing effects are clearly 
reflected in the staggering behavior of A versus A as seen on the right column of 
Fig. 3. In the empirical formulae, A is close to zero for odd-even nuclei, positive 
for even-even nuclei and negative for odd-odd nuclei. However, we see that both 



P(E X ) 




(11) 



7 6 
> 



1 1 1 1 
Mn 


i i i i i i i 
X 

></■ 


i i i i 

. Fe 


i i i i ^ i i 


: ^ 




1 1 1 1 
- Co 


—\ — i— i — i — i — 

x 

,m / 
* / 







50 52 54 56 58 60 

A 



- 2 A / ^' V 

-4 ' 




- < -2 - 



-4 





-2 - 
-4 



Mn 



X 



1 1 1 1 — ± — I 1 1 1 1 h 



I I I I I I I I I I I 



Co 



i \ 



x* 



50 52 54 56 58 60 

A 



FIGURE 3. The single-particle level density parameter a (left column) and the backshift pa- 
rameter A (right column) for Mn, Fe and Co isotopes. The solid squares are the SMMC values 
(obtained by fitting the SMMC level densities to the BBF (11)). The x's are the experimental 
results taken from the compilations of [4] (assuming rigid body moment of inertia), except for 
58 Co and 59 Co where we used the values quoted in [5]. The solid lines are the empirical formulae 
of [14]. Taken from Rcf. [12]. 



the experimental and SMMC values of A can differ significantly from zero for the 
odd-even nuclei. The SMMC values of a and particularly A are generally in better 
agreement with the experimental results than the empirical values. For some of the 
odd-odd manganese isotopes we observe discrepancies between the SMMC values 
of a and the experimental data. However, the lower values of a for these manganese 
isotopes are compensated by corresponding lower values of A. Consequently, the 
discrepancies in the level densities themselves are less significant for E x < 10 MeV. 

To demonstrate the T z - dependence of level densities we show in Fig. 4 the level 
densities of two odd-odd A = 54 (Mn and Co) and three odd-even A = 55 nuclei 
(Mn, Fe and Co). The empirical formulae predict similar level densities for the two 
odd-odd nuclei as well as for the three odd-A nuclei: the values of a are similar if 
the mass A is the same; A ~ for odd-A nuclei; and A are approximately the same 
for odd-odd nuclei. However the SMMC level densities of these nuclei (symbols) 
are seen to be quite different from each other. We also see that the experimental 
level densities (dashed lines) are in good agreement with the SMMC densities. 

In conclusion, we have described a particle-number reprojection method in the 
shell model Monte Carlo method. With this reprojection technique we can calculate 
the thermal properties for a series of nuclei using Monte Carlo sampling for a single 
nucleus. Level densities of odd-A and odd-odd nuclei are calculated despite a new 
sign problem introduced by the projection on an odd number of particles. 
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FIGURE 4. The level densities of two odd-odd A = 54 nuclei: 54 Mn and 54 Co (left), and three 
odd-even A = 55 nuclei: 55 Mn, 55 Fe and 57 Co (right). The symbols are the SMMC level densities 
and the dashed lines are the experimental level densities. 
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